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Abstract 

We consider a special case of the three dimensional Vlasov-Poisson 
system where the particles are restricted to a plane, a situation that is 
used in astrophysics to model extremely flattened galaxies. We prove 
the existence of steady states of this system. They are obtained as 
minimizers of an energy- Casimir functional from which fact a certain 
dynamical stability property is deduced. Prom a mathematics point 
of view these steady states provide examples of partially singular so- 
lutions of the three dimensional Vlasov-Poisson system. 

1 Introduction 

In astrophysics the time evolution of large stellar systems such as galaxies is 
often modeled by the Vlasov-Poisson system: 

dtf + v-VJ-V.U-VJ = 0, 
AU = A7cp, lim U{t,x) = 0, 



p{t,x) = / f{t,x,v)dv. 
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Here f = f(t,x,v)>0 denotes the density of the stars in phase space, tGlR 
denotes time, x,f gH^ denote position and velocity respectively, p is the 
spatial mass density, and U the gravitational potential. The only interaction 
between the stars is via the gravitational field which the stars create collec- 
tively, in particular, collisions are neglected. When modeling an extremely 
fiattened galaxy the stars can be taken to be concentrated in a plane (the 
(xi,a;2)-plane). The corresponding potential which is given by the usual inte- 
gral representation induces a force field which accelerates the particles only 
parallelly to the plane, and the Vlasov-Poisson system takes the form 



dtf + v-VJ-V,U-\/J = 0, telR, x,veM', (1.1) 

U{t,x) = -f ^^dy, (1.2) 
jir2 \x — y\ 



pit^x)= / Jit,x,v)dv. (1.3) 

Note that from this point on, x,f GlR^. The three dimensional phase space 
and spatial densities are given as 

f{t,X,X3,V,V3)=f{t,X,v)S{x3)5{v3) 

and 

p{t,X,X3)=p{t,x)5{x3) 

where S denotes the Dirac distribution. It should be emphazised that the sys- 
tem (|1.1|) , (|1.2|) , (|1.3|) is not a two dimensional version of the Vlasov-Poisson 
system but a special case of the three dimensional system with partially sin- 
gular phase space density. In the present paper we are concerned with the 
existence of steady states of this system and with their stability properties. 
There are a number of aspects which make this problem interesting. Al- 
though such fiat solutions of the Vlasov-Poisson system occur as models in 
the astrophysics literature, cf. @, |^, we know of no mathematical investiga- 
tion of this situation. The fact that the distribution function is singular in 
the a:3-direction, or, alternatively, that the two dimensional Vlasov equation 
is coupled to a potential with the three dimensional l/|a:|-singularity, makes 
this problem mathematically nontrivial. We refer to [0, where solutions of 



the Vlasov-Poisson system which are measures are treated in the one dimen- 
sional case; an extension of these results to higher dimensions is not known. 
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Finally, the method that we employ to study the existence and the stabil- 
ity properties of steady states was recently used in a spherically symmetric, 
regular, three dimensional situation in |TT| . The present paper demonstrates 



that this method extends beyond the case of spherical symmetry, although 
this assumption played an important role in |1TT| . 

To see how steady states of the system ( |1.1| ), ( |1.2| ), ( |1.3| ) can be obtained, 
note first that if U = U{x) is time independent, the particle energy 

E = ^\v\^ + U{x) (1.4) 

is conserved along characteristics of (|1.1|) . Thus any function of the form 

f{x,v) = <PiE) (1.5) 

satisfies the Vlasov equation. We construct steady states as minimizers of 
an appropriately defined energy- Casimir functional. Given a function Q = 
<5(/)>0, />0, we define 

V{f):= JlQif)dvdx + ^J J \v\'fdvdx + ypfUfdx. 

Here / = f{x,v) is taken from some appropriate set J-'m of functions which in 
particular have total mass equal to a prescribed constant M, pf denotes the 
spatial density induced by / via ( |1.3|) , and Uj denotes the potential induced 
by Pf via ( |1.2| ). If one can show that the functional V has a minimizer, 
then the corresponding Euler-Lagrange equation turns out to be of the form 
(11.51). An alternative way to obtain steady states would be to substitute ( |1.5| ) 
into ( |1.3| ) so that p would become a functional of U, and it would remain to 
solve ( p.. 21) which becomes a nonlinear integral equation for U. This route is 
followed for example in |l|] for the regular three dimensional problem. The 
major difficulty then is to show that the resulting steady state has finite mass 
and compact support — properties which are essential for a steady state to 
qualify as a physically viable model — , and this problem has been dealt with 
for the polytropic ansatz 



f{x,v) = {E-E^ 







where Eo is a constant, — l</i<7/2, and (■)+ denotes the positive part. Our 
approach has the advantage that finiteness of the total mass and compact 
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support are built in or appear naturally, and these properties do not depend 
on a specific ansatz like the polytropic one. Furthermore, the fact that the 
steady state is obtained as a minimizer of the functional T) implies a certain 
nonlinear stability property of that steady state. 

The paper proceeds as follows. In the next section the assumptions on 
the function Q which determines our energy-Casimir functional are stated, 
and some preliminary results, in particular a lower bound of V on Tm, are 
established. The main difficulties in finding a minimizer of T> arise from 
the fact that V is neither positive definite nor convex, and from the lack 
of compactness: Along a minimizing sequence some mass might escape to 
infinity. However, using the scaling properties of V and a certain splitting 
estimate we show that along a minimizing sequence the total mass has to 
concentrate in a disc of a certain radius Rmi depending only on M. In [|ll| the 
corresponding argument required the assumption of spherical symmetry. In 
the present paper we only require axial symmetry with respect to the xs-axis. 
The corresponding estimates are proved in Section 3 and are used in Section 4 
to show the existence of a minimizer. It is then straight forward to show that 
the Euler-Lagrange equation is (equivalent to) ( |1.5|) , thereby completing the 
existence proof for the steady states. The resulting stability property of such 
steady states is discussed in Section 5. Since we have to restrict our functions 
in the set J^m to axially symmetric ones stability holds only with respect to 
such perturbations, and also perturbations transversal to the (xi,X2)-plane 
are not covered. Moreover, the stability result is only conditional in the 
sense that so far no existence theory for the initial value problem for the flat 
Vlasov-Poisson system is available. To obtain a complete stability result, 
global existence of solutions which preserve the energy-Casimir functional P 
would be needed, at least for data close to the steady states. In the last 
section we briefly discuss the regularity properties of the obtained steady 
states. 

We conclude this introduction with some references to the literature. In 
the regular three dimensional situation the existence of global classical solu- 
tions to the corresponding initial value problem has been shown in |]T^, cf. 
also 0, [15|, ^ . The existence of steady states for the case of the polytropic 
ansatz was investigated in ^ and [Q. We refer to @] for contributions to 
the stability problem in the astrophysics literature. As to mathematically 
rigorous results on the stability problem, we mention [1T0| , p!T|| for applications 
of the present approach in the regular, three dimensional case, cf. also 
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An investigation of linearized stability is given in 0. For the plasma physics 
case, where the sign in the Poisson equation is reversed, the stability problem 
is much easier and better understood. We refer to |T^, [T^, [T^. A plasma 
physics situation with magnetic field is investigated in 0. 

2 Preliminaries; a lower bound for V 

We first state the assumptions on Q which we need in the following: 

Assumptions on Q: For Q eC^{[0,oo[), Q>0, and constants Ci, . . . ,6*4 > 0, 
Fq > 0, and < /ii, /i2, /U3 < 1 consider the following assumptions: 

(Qi) g(/)>Ci/i+i/''\ />Fo, 

(Q2) g(/)<C2/i+i/'-, 0</<Fo, 

(Q3) Q{Xf)>X'+y^-^Qif), />0, 0<A<1, 

(Q4) Q"(/)>0, />0, andg'(0) = 0. 

(Q5) C3Q"{f) < Q"{Xf) < CiQ'\f) for / > and A in some neighborhood of 
1. 

The above assumptions imply that Q' is strictly increasing with range [0,oo[, 
and we denote its inverse by g, i. e., 

g'(g(e)) = e, e>0; (2.1) 

we extend q by g(e) = 0, e < 0. 

Remark: The steady states obtained later will be of the form 

Ux,v) = q{Eo-E) 

with some < and E as defined in (|T1]). If we take Q{f) = /^+^/^, / > 0, 
this leads to the polytropic ansatz, and such a Q satisfies the assumptions 
above if < yU < 1. If we take 

Q(/)=Ci/^+^/^^+C2/^+^/^^ (2.2) 
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with 0</ii,/i2<l and constants Ci,C2>0 then again the above assump- 
tions hold, but q is not of polytropic form. Due to the assumption of axial 
symmetry which we will have to make for other reasons, 

Ls = XiV2-X2Vi, 

the Xs-component of angular momentum, is conserved along characteristics 
as well. It would be a purely technical matter to allow for depence on L3 of a 
type where for example the constants in (|2.2| ) could be replaced by functions 
of L3 which are bounded and bounded away from 0. We refer to |jll|] for the 
necessary modifications. 

For a measurable function f = f{x,v) we define 

Pfi^) '■= j fix,v)dv 
and ^ 

as to the existence of this convolution see Lemma § below. Then define 

1 rr. ,2. 



Ekinif) '=2]] 1^1 fix,v)dvdx, 
EUf) ■■= \ j PA^)Uf{x)dx = -\j j PM^dxdy, 

C(/):= JjQ{f{x,v))dvdx, 

V{f):=EUf)+C{f), 
V{f):=V{f)+E^M)- 

The sum £^kin(/) + -Epot(/) is the total energy corresponding to /, a conserved 
quantity for the time dependent problem, and the same is true for C, a 
Casimir functional which corresponds to the conservation of phase space 
volume. T> is the energy-Casimir functional, and V is the positive part of 
that functional. We will also use the notation Up and -E'pot(p) if P = p{x) is 
not necessarily induced by some f = f{x,v). The following two sets will serve 
as domains of definition for the energy-Casimir functional V: 

^m:={/gI^'(]R')|/>0, j Jfdvdx = M, P(/)<oo|, (2.3) 
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where M > is prescribed, and 

'■— {/ ^ I / is axially symmetricj. (2.4) 

Here axial symmetry means that 

f(Ax,Av)^f{x,v), x,veWi^, AeS0{2). 

When viewed as a function on the effective phase space IR^, / is spherically 
symmetric, but when viewed as a function over the full phase space H^, / is 
only axially symmetric. The induced potential does not share the properties 
of spherically symmetric potentials which is why we prefer the phrase axially 
symmetric. We do not restrict ourselves to the set JF^ from the beginning 
in order to point out where exactly the symmetry is needed. 

The aim of the present section is to establish a lower bound for D of a 
form that will imply the boundedness of V along any minimizing sequence. 

Lemma 1 Let (Ql) hold and let ni = l + fii. Then there exists a constant 
C > such that for all f G J-'m, 

j pY^''''dx<C{M + V{f)). 

Proof : We split the v integral into small and large v's and optimize to obtain 
the estimate 

//■,./ \ 2mi/(4+2m) / . n2/(4+2mi) 

Pf(x)<ci^J f'+'/'^^dvj [j\v\'fdvj 
By definition of ni and assumption (Ql) we find 

J pY^''''dx <c(^J f^+^/^''dvdx + J \v\^ fdvdx+ 

<c(^fI^^'^' j fdvdx+^ jQ{f)dvdx + J \v\^ fdvdx+ 

and by definition of V this is the assertion. □ 
Note that 3/2 < 1 + 1/ni < 2, and since by definition pfE L^(IR^) for / e 
J^M, we have p/ e L'^/^i'iR^) for feJ^M- 
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Lemma 2 //p G L^^^{JR'^) then Up E L^(IR^), and there exists a constant C > 
such that for all p e L^/^{Wi^) , 

\\Up\\,<C\\p\\,/„ -i?pot(p)<C||p||J/3. 

Proof : The assertion follows from generalized Young's inequality p. 32], 
since 1/| • | gL^(IR^), the weak L^-space, and from Holder's inequality. □ 

Combining the previous two lemmata yields the desired lower bound of 
V over the set J^m'- 

Lemma 3 Let Q satisfy assumption (Ql). Then 

PM: = inf{P(/)|/e^M}>-oo, 
and there exists a constant Cm>0 depending on M such that 

and for any minimizing sequence (fn) CJ-'m ofD we have 

V{fn)<CM, nelN. 

Proof : If we interpolate the L^/'^-norm between the L^-norm and the 
norm and apply Lemma |] we find 

/ pfdx < Cm (/ < Cm {l+V{f)r'' . 

Thus by Lemma ^ 

2— - 

V{f)>V{f)-CM{l + V{f)) ^'>P(/)-Cm(i + P(/)"^/2). 

Since ni<2 the rest of the lemma is obvious after possibly choosing Cm 
larger. □ 

In later sections we will have to assume axial symmetry, and we will need 
the fact that 

Pf,:=inf{P(/)|/G^f,}>-oo, 
which of course follows from the previous lemma. 
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3 Scaling and Splitting 



The behaviour of T> and M under scahng transformations can be used to 
relate the Pm's for different values of M: 

Lemma 4 Let Q satisfy the assumptions (Q1)-(Q3). Then —oo<Vm<0 
for each M>0, and for all < Mi < M2, 

where a = 1/(1 — /is) > 0. The same assertions holdforVfj instead of Dm- 

Proof : Given any function f{x,v), we define a rescaled function f{x,v) = 
af{bx,cv), where a, b, oO. Then 

J J fdvdx = ab~'^c~'^ J J fdvdx (3.1) 

and 

V{f) = b-^c'^Ciaf) + ab-^c-^E^Uf)+a'b-'c-^Ep,tif). (3.2) 

Proof ofVM<0: Fix some /GJFf with compact support and f<Fo, and 
let a = Mb'^(? so that f ETfj. The last term in P(/) is negative and of the 
order 6, and we want to make this term dominate the others as 6 ^ 0. Choose 
c = b~^l'^ so that a = Mb^~^ ^ and assume that a<l so that af<Fo. By (Q2), 

'Dif) < C (^6(2-7)^2 ^ 57^ _cb 

where C, C >0 depend on /. Since we want the last term to dominate as 
b-^0, we need 7>1 and (2 — 7)//i2>l, and, in order that a<l as b^O, 
also 7 < 2. Such a choice of 7 is possible since fi2 < 1, and thus V^f) < for 
b sufficiently small. 

Proof of the scaling inequality: Assume that / G J-'m2 and / G J^Mi so that by 

ab-^c-' = ^=:m<l. (3.3) 
M2 - ^ ^ 

By dU) and (Q3), 

Vif)>ma'/'^'Cif) + mc-^EUf)+m'bEpM) 
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provided a < 1 . Now we require that 

Together with ( ^.31 ) this determines a, b, c in terms of m. In particular, 
as required and 

P(/)>m'+^P(/). 

Since for any given choice of a, b, c the mapping / 1— >■ / is one-to-one and onto 
between J^M2 and J-'mi as well as between J^^^ and J-'§^_^ the scaling inequality 
follows. □ 
The following two lemmata are crucial in proving that along a minimizing 
sequence the mass concentrates in a certain ball. It is here that we need the 
additional symmetry assumption and where the estimates become more in- 
volved than in the regular spherically symmetric case. The aim is to estimate 
the effect on T> of splitting the matter distribution into a part inside a ball 
Bji of (large) radius R about and a part outside. 

Lemma 5 There exists a constant C > such that for every p&L^n 
L'^^^(]R'^) which is nonnegative and axially symmetric, i. e., p{x)=p{\x\), 
and every R>0 the following estimate holds: 

-f p(x)Up(x)dx<CR-^/'^\\ph,.. f p(x)dx. 

J\x\>R ' J\x\>R 

Proof : Due to the symmetry of p the potential is given by 

U,{x) = U,ir) = -4 r ^pis) K f ^) ds 
Jo r + s \r + s J 

where the elliptic integral K is defined as 

dd) /-i dt 
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We need to estimate the singularity in K: 



Jo 



dt 1 , 1 + 



<C(l-ln(l-0), 0<e<l. 
Substituting for ^ yields 

for < r < s; the case r > s is analogous. Thus 

^/V^\ -ln(l-[r,s])), r,s>0 (3.4) 

\r + s J 

where 

( T S 

[r, s] : = min 



s r 



Now 



where in Ji the variable r ranges in [0,2s] and in I2 it ranges in [2s,oo[. Using 
and Holder's inequality we find 

h<C p{s) rp{r){l-\n{l-[r,s\))drds 



R 



<C\\p\\,/,l^pis) (^'%(l-ln(l-[r,s]))V)'^'rfs 
<C\\ph/J s'/'p{s)ds<C\\p\\,/,R-'/' sp{s)ds- 
note that with a = r/s, 

J\ (1 - \n{l - [r, s])y dr = J\(l-\n{l -a)y da 

+ a (r)y da 

= Cs\ s>0. 
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The second term is much easier to estimate: For r>2s we have — ln(l — 
s/r) <ln2, and by Holder's inequahty, 

/■oo /"oo ^ poo poo 

h ^ C sp(s) / p{r)drds<C / sp{s)ds / p{r)dr 

JR J2s r + s Jr Jr 

poo 

<C\\p\\^/.^R-^/^ sp{s)ds. 

J R 

Together with the estimate for Ji this completes the proof. □ 
Lemma 6 Let Q satisfy the assumptions (Q1)-(Q3) and let f eJ-'^- Then 

W)-^^ff>f^/ ffdvdx-^]f ffdvdx,R>0, 

V ^ J\x\<RJ y/R)J\x\>RJ 

where the constant < depends on a from Lemma ^ and Cm > depends 
on M. 

Proof : Let Br denote the ball of radius R about in IR^, let Ibrxr^ be the 
characteristic function of Bji x IR^, 

/i = 1_BbxR2/, /2 = / — /l, 

and let pi and Ui denote the induced spatial densities and potentials respec- 
tively, i = l,2. We abbreviate X = Jf2- Then 

T^{f)=V{fi)+V{f2) + \jUipidx + ]^j U2P2dx + J Uip2dx 

>vI,_^+vI-CmR-^'^\ 

since /i G J^m_x and /2 G J-'^. To estimate the "mixed term" in the potential 
energy we have used Lemma ^; note that for / G J-'m, 1 1 P/ 1 1 4/3 is bounded by a 
constant depending only on M, cf. Lemma |l|. Since a > 0, there is a constant 
Ca<0, such that 

(l-x)^+" + x^+"-l<C„(l-a;)x, 0<a;<l. 
Using Lemma ^ and noticing that Vf^ < we find that 

V{f)-Vlj > [(l-A/M)i+" + (A/M)i+"-l]pf^-CMi?"^/^A 
> C^VI,{1-X/M)X/M-CmR'^^^X 
which is the assertion. □ 
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4 Minimizers of T> 

Before we show the existence of a minimizer of T> over the set JF^^ we use 
Lemma ^ to show that along a minimizing sequence the mass has to concen- 
trate in a certain baU: 



Lemma 7 Let Q satisfy the assumptions (Q1)-(Q3), and define 

Rm 



(2MCmV 



where Cq < and Ca/ > are as in Lemma // (/„,) C Tfi is a minimizing 
sequence ofV, then for any R> Rm, 



hm / / fndvdx = 0. 
'^-'°°J\x\>rJ 



Proof : If not, there exist some _R> Rm, A > 0, and a subsequence, called (/„) 
again, such that 

lim / fndvdx = \. 



l\x\ 

For every ri G IN we can now choose Rn> R such that 



:= / / fndvdx= \ I I fndvdx. 

J\x\>RriJ IJ\x\>rJ 



Then 



lim / / fndvdx = lim A„ = A/2>0. 
Applying Lemma ^ to Bfi^ we get 



M2 '^-^ '-^ V M2 '^-^ ^ 



Cm \ 


1 Xn> 




VR^ / 


[ M2 






^ CaT^M 


Vrj 


, 2M 



/ CaVli Cm \ ^ ( Cal^ii Cm \ X 

^[-m^^^-^^^^-7r)2-[-^~7r)2 

as — > oo, since < A < M. By definition of Rm the expression in the paren- 
thesis is positive for R>Rm, and this contradicts the fact that (/„) is a 
minimizing sequence. □ 

As a further prerequisite for the existence proof of a minimizer we estab- 
lish a compactness property of the potential energy functional: 
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Lemma 8 Let (p„) cL^/^nL^(IR^) be bounded and axially symmetric with 
Pn^Po weakly in L'^/^(IR^), n^oo. 

Then 



EpotiPn- Po)^0 and EpotiPn) ^ Epotipo), n^oo. 

Proof : We consider the convergence of -Epot(Pn — Po) first. By Lemma |^ and 
Holder's inequality 



\Epot{Pn- Po) \ < 



iUp„-Up,,){pn-po)dx 



< l|Pn-po|l4/3l|f^p„,i?-f^po,i?llL4(B^) + -^, 



for any R>0, where 



Up,r{x) :=- I f^^^ dy, x e IR^. 



f\y\<R \x-y\ 

Thus it suffices to show that for R>0 fixed the mapping 

T:L'/\Bn)3p^lB,{p*k)eL\JR') 

is compact where k := 1^2^? 1/1 ' |j ^ote that we may cut off the Green's func- 
tion as indicated since only x,y with |x|,|?/|<i? need to be considered. 
We use the Frechet-Kolmogorov criterion to show that T is compact. Let 
S <zL^^'^{Br) be bounded. Then TS is bounded in L'^{W{?) by Lemma ^ 
Since the elements in TS have a uniformly compact support it remains to 
show that for p G 5, 

{Tp)h^Tp in L^(1R2), /i^O, 
where gh:= g{- + h), HeIR^. But by Young's inequality, 

\\{Tp)h-Tp\\^<\\p*{kh-k)\\^<\\p\\yJkh-k\\^y^^O 
uniformly on S as h^O, since kEL^'^^'^iJR^). Since 

^pot (Pn) - Epat (Po) = - y f^po (pn - Po) " ^pot {pn " Po) 

the fact that Up^ gL^(IR^) together with the weak convergence of p„ implies 
the remaining assertion. □ 
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Theorem 1 Let Q satisfy the assumptions ( Ql )-( Q4 ), and let (fn) C J^^j 
be a minimizing sequence ofT>. Then there is a minimizer /oGJF^j and 
a subsequence {fn,,) such that V^fo) = Vf.j, supp /o C -Bj?^^ x IR^ with Rm as 
in Lemma ^ and fn^^fo weakly in L^~^^^^^ (JR'^) . Furthermore, i?pot(/rafc~ 
/o)-0. 

Proof : By Lemma 0, (V{fnj) is bounded. Let pi = l + l//ii. Then the 
sequence (/„) is bounded in Lpi(IR^) by assumption (Ql). Thus there exists 
a weakly convergent subsequence, denoted by (/„) again, i. e., 

/„-/o weakly in Lf^(lR^). 

Clearly, /o>0 a. e., and /q is axially symmetric. Since by Lemma ^ 

M = lim / / fndvdx+ lim / / fndvdx 

J|x|<_Ri J\v\<R2 "'-*°°J\x\<B^i J\v\>R2 

r r C 
< lim / / fndvdx + —^ 

" \x\<RiJ\v\<R2 K2 



n^oo 



where Ri> Rm and -R2 > are arbitrary, it follows that 

/ I fndvdx = M 

J\x\<RiJ 

for every Ri > Rm- This proves the assertion on supp/o and ///o = M. Also 
by weak convergence 

1 1 < liminf J J \v\'^ fndvdx < 00. (4.1) 

By Lemma (p„) = (p/„) is bounded in 0^^) where ni= ni + 1. After 

extracting a further subsequence, we thus have that 

Pn^Po:=P/o weakly in L^/^(1R^), 

and Lemma ^ implies the convergence of the potential energy term. 

It remains to show that /o is actually a minimizer, in particular, 
V{fo)<oo so that foEj-'lj. By Mazur's Lemma there exists a sequence 
((7n) C LP^(IR^) such that gn-^fo strongly in Lp^(IR^) and gn is a convex 
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combination of {fk \ k>n}. In particular, fo a- e. on IR^ By (Q4) the 
functional C is convex. Combining this with Fatou's Lemma implies that 

C(/o) <liminfC(f7„) <limsupC(/„). 

Together with ( [4.1|) this imphes that 

P(/o)<JimP(/„)<oo; 

note that lim„^oo'P(/n) exists. Therefore, 

V{fo)=V{fo)+E^Mo)<}im{V{U)+E^Mn)) = Vi, 
and the proof is complete. □ 

Theorem 2 Let Q satisfy the assumptions (Q1)-(Q5), and let JqEJ-'m be a 
minimizer ofD. Then 

fo{x,v) = q{EQ — E) a. e. on IR"^, 

where 

2' 



E = l\v\^ + Uo{x) 



Eo = j^ll{Q'{fo)+E)fodvdx<0, 



Uq is the potential induced by /q, and q is as defined in ( \2.1\) . 



Note that UQ = Uf^ by construction, and /o is a function of the particle energy 
only and thus a steady state of the system (|l.lj ), ( |1.2|) , ( |1.3|) . The regularity 



of f/o and thus the sense in which /o satisfies the Vlasov equation ( p..l[ ) is 
investigated in the last section. 

Proof : Let /o be a minimizer. For fixed e > let rj: IR^ — > IR be measurable, 
with compact support, axially symmetric, and such that 

|?7| < 1, a. e. on IR"^, r/ > a. e. on IR'^ysupp/o, 

and ^ 

e < /o < - a. e. on supp/o flsuppr;. 
e 
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Below we will occasionally argue pointwise on IR'^ so we choose a represen- 
tative of /o satisfying the previous estimate pointwise on supp/o flsupp?]. 
For 

0<h<— — 
2(l + ||r/||i) 

we define 

gih) = M, + 



ll^^ + /o||i' 

This defines a variation of /o with g{h) e and 5'(0) = /o; note that 

M-U\\hrj + fo\\i<M + ^. 
We expand V{g{h)) —'D{fQ) in powers of h: 

V{g{h))-V{fo) {Q{g{h))-Q{fo))dvdx+y J \v\'{g{h) - fo)dvdx 

+ 11 Uo{g{h)-U)dvdx + E^,,{g{h)-U). (4.2) 

Since g{h) >0 on IR^, g{h) is differentiable with respect to /i, and we write 
g'{h) for this derivative. Note that both g{h) and g'{h) are functions of 
{x,v) elR^, but we suppress this dependence and obtain 



Now 
and 

so that on IR^, 



9'{^)-V-^[j jvdvdx^U (4.3) 

\9"{h)\<C{\n\ + h) 
\g{K)-fo-hg'm<ChWr^\+fo)-, 
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in the following, constants denoted by C may depend on /o, ?7, and e but 



never on h. We can now estimate the last three terms in ( 4.2 



j J \v\^{g{h)-fo)dvdx = h J J \v\^g'{0)dvdx + O{h^ 

J JUo{g{h)-fo)dvdx = hJ J Uog'{0)dvdx + O{h^) 
\EMgih)-fo)\<C\\pg^H)-Po\\l,<Ch^ 

For the last estimate we used Lemma |^ and the fact that 
\Pgih){x)-po{x)\<Ch / {\r]\+fo){x,v)dv. 



(4.4) 

(4.5) 
(4.6) 



It remains to estimate the first term in (4^). Consider first a point {x,v) G 
supp/o with /o(x,f)>0. Then 

Q{9{h)) - Q(/o) = Q'{fo){9{h) - /o) + ^g"(r) {g{h) - Uf 

= hQ\fo)g'iO) + h'^Q'ifo)g"ie) + \Q"ir)igih) - Uf 

where r lies between g{h) and /o and 9 lies between and h\ both r and 9 
depend on (x,w). Thus 

Q(^W)-Q(/o)-/ig'(/o)^'(0)|<Cg'(/o)(|r/|+/o)/^' + CQ"(r)(|r/p + /o2)/.2. 
On supp/o we have 

\h<g{h)<2h 

provided 0<e<eo with eo>0 sufficiently small. Thus r lies between /o/4 
and 2/o, and by iterating (Q5) a finite, /^-independent number of times we 
find 

Q"{r)<CQ"{h). 

By (Q3) and (Q5), 



and thus 



Q(^W)-Q(/o)-/iQ'(/o)^?'(0) <CQ{h)h^ + C\7^\h 



,2. 



here we used the continuity of Q' and Q" and the fact that e < /o < 1/e on 
suppr^risupp/o. The above estimate holds for any point {x^v) Gsupp/o with 
fo{x,v) >0. Now consider a point {x,v) with /o(a;,f) =0. Then 

9{h)=M <C\v\h 
\\hv + fo\\i 

so that by (Q4) and (Q2), 

Q{9{h))-Q{fo)-hQ'{h)9'{0)\ = Q{g{h)) < Q{Ch\r^\) 

for h>0 sufficiently small. Thus 

Qi9ih))-Qifo)-hQ'ifo)g'iO) dvdx<Ch'^' (4.7) 



for some 6>0. Combining ( ^ , (U), (^), and ( ^ ) with the fact that /o 
is a minimizer we find 

^<V{g{h))-V{h) = hj j {Q\f,) + ]^\v\^ + U,^g'{Q)dvdx + 0{h^+') 



for all h>{] sufficiently small. Recalling ( ^.3| ) and the definitions of E and 
£^0 this implies that 



/ J{Q'Uo)+E-Eo)vdvdx>0. 



Recalling the class of admissable test functions rj and the fact that e > is 
arbitrary, provided it is sufficiently small, we conclude that 

E — Eo>0 a. e. on IR'' \ supp /o 

and 

Q\fo)+E-Eo = a. e. on supp /o. 
By definition of q — cf. ( p.l| ) — this implies that 

fQ{x,v) =q{EQ — E) a. e. on IR^. 

Since po has compact support and lim2,^oof^o(a;) = we conclude that Eq < 0. 

□ 
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5 Dynamical Stability 



We now discuss the dynamical stability of /q. As noted in the introduction 
the existence of solutions to the initial value problem for the system (|1 . 1|) , 
L2D , (|1.3| ) is open. In the following we therefore have to assume that for 



initial data in some (reasonably large) set X C JF^^ the system has a solution 
f{t) with f{t)eTli and V{f{t))=V{f{0)), t>0; classical solutions of the 
regular three dimensional Vlasov-Poisson system have these properties. The 
considerations below are only formal, and we emphasize this fact by not 
stating any theorems but only giving the stability estimates. First we note 
that for / G J-'m, 

V{f)-V{U)=d{fJo)+E^M-fo) (5.1) 

where 



d{fJo)= Q{f)-Q{fo) + {E-Eo){f-fo] 



dvdx. 

Next we observe that d{f ,/())> 0, JeJ-'m- For E — Eq>0 we have /o = 0, 
and thus 

W) - Wo) + (i^ - i^o) (/ - /o) > g(/) > 0. 

For E-Eo<0, 

Qif) - QUo) + {E- E,) if - /o) = ^g"(/) (/ - > (5.2) 

provided / > 0; here / is between / and /q. If / = 0, the left hand side is still 
nonnegative by continuity. Now let Q satisfy the assumptions (Q1)-(Q5) and 
assume that the minimizer /o is unique in JF-^. Then we obtain the following 
stability estimate: 

For every e > there is 6>0 such that for any solution f{t) of the fiat Vlasov- 
Poisson system with /(O) G X , 

rf(/(0),/o) + |i^pot(/(0)-/o)|<5 

implies 

difit),fo) + \E^otifit)-fo)\<e, t>0. 

Assume this assertion were false. Then there exist eo > 0, t„ > 0, and /„(0) G 
X such that 

C^(/n(0),/o) + |i^^pot(/n(0)-/o)| = ^ 
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but 

diUQJo) + |^pot(/n(tn) - /o)| > eo > 0. 

From ( |5.1| ), we have lim„_»oo2^(/n(0)) =2^ff. Since V^f) is invariant under 
the assumed Vlasov-Poisson flow, 

hm P(/„(t„))= hm V{UO))=Vl,. 

Thus, {fn{tn)) (^^ii is a minimizing sequence of V, and by Theorem |1| , we 
deduce that — up to a subsequence — Epot{fn(tn) — fo) Again by (|5.1|) , 
d{fn(tn),fo)^0, a contradiction. 

Provided the assumed global Vlasov-Poisson flow is such that in addition 
II /Wlloo = 11/(0) Hoc. ^>0' and that Q is such that 

Ci:=inf{Q"(/)|0</<Co}>0 

for some constant Cq > ||/o|loo; then for /(0)<Co one obtains the stability 
estimate 



//, Q{f{t))dvdx + ^ 11 \f(t)-fo\'dvdx+\E,Mit)-fo)\<e. 



/R^Xsupp/o 2 J Jsupp/o 

This follows by estimating Q" in the expansion (|5.2|) from below. 

If the minimizer /o of D is not unique (and not isolated) in jF|j, then a 
solution starting close to /o — in the sense of the above measurement for the 
deviation — remains close to the set of all minimizers in Tf^^. In the regular, 
three dimensional case uniqueness of the minimizer can be shown for the 
polytropic ansatz, cf. |TT . 



6 Regularity 

So far the steady states obtained in Section 4 satisfy the Vlasov-Poisson 
system ( |lT1) , (|1.2|) , ( p.. 3D in a rather weak sense, in particular, the potential 
need not be sufficiently regular for characteristics of the Vlasov equation to 
exist so that the precise meaning of /o being a function of an invariant of 
the particle trajectories is questionable. The present section will remedy this 
under some very mild additional assumptions: 
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Theorem 3 Assume that Q satisfies conditions (Q1)-(Q5), and in addition 

Q'U)>Cif''\ f>Fo. 

Let fo^^ij ^6 ^ minimizer ofD as obtained in Theorem 1, and po, Uq the 
induced spatial density and potential respectively. Then po, Uq G C^(JR^), and 
the first derivatives of Uq are Holder continuous. If the function q defined in 
(\2.1\ ) is locally Holder continuous, then UoEC^iJR^), and the second deriva- 
tives of Uq are Holder continuous. 

Proof : As a first step we wish to show that Uq and po are bounded. RecaU 
that 

-Uo{r)=4 ^p,{s)K{Ods = h + l2,r>0, 
Jo r + s 

where in Ii the variable s ranges in [0,2r] and in I2 it ranges in [2r,oo[. Using 
the estimate (p.4|) and the fact that po E L^^^iJR"^) we find 

C /-Sr 



h< - f ^ spo{s)(l-\n{l-[r,s]))ds 
r Jo ^ ' 

< 711^0113/2(7^ rs(l-ln(l-[r,s])) dsj <Cr-^l\ 



For s > 2r the elhptic integral K{^) is bounded, and again by Holder's in- 
equality we immediately obtain the same estimate for I2 so that 

|f/o(r)|<Cr-l/^ r>0. 



Next we know that 

, , ( { ^ 1 . s\ , I 27r / 

PoKr) 



27T q{Eo-E)dE ,Uo{r)<E, 



0, 

,f/o(r)>^o. 



The additional assumption on Q' implies that there are constants C > 0, eo > 
such that 

q{t)<Ce^\ e>eo, 

and this implies that 

Po(r)<Cr-(^l+^)/^ r>0. 
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Since we know that po has compact support, it follows that po G L^(IR^). We 
may now repeat the estimate for Uq and obtain Ii < Cr^l"^ and /2 < C so that 
?7o and thus also po are bounded. 

For the rest of our argument we rely on the regularity properties of po- 
tentials generated by single layers. Firstly, the boundedness of po implies 
that [/q is Holder continuous, cf. page 42]. The relation between po and 

immediately implies that po shares this property. This implies that 
has Holder continuous first order derivatives, a fact known as Ljapunov's 
Theorem, cf. |0, pages 66, 67]. Since 

p'o(r) = -27rg(Eo-t/o(r))[/^(r) 

Po is continuously differentiable. 

If q is locally Holder continuous, then po will have Holder continuous 
first order derivatives; note that Eo — UQ{r) ranges in a bounded interval for 
rG[0,oo[ so the local Holder continuity of q suffices. We can now apply 
Ljapunov's Theorem again and obtain the remaining assertions. □ 

We remark that above we considered Uq as a function on the {xi,X2) plane. 
Of course the definition ( |1.2| ) makes perfect sense on all of IR^, and as long 
as we consider only derivatives parallel to the (xi,X2) plane all the regularity 
assertions for Uq hold on the whole space IR^. However, it is well known that 
the derivative of Uo perpendicular to the plane has a jump discontinuity on 
the plane. 
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